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(N ■ Abstract 

We present the analytic Lifshitz solutions for a scalar field model non minimally coupled 
^ . with the abelian gauge field in N dimensions. We also consider the presence of cosmological 

constant A. The Lifshitz parameter z appearing in the solution plays the role of the Lorentz 
breaking parameter of the model. We investigate the thermodynamical properties of the 
solutions and discuss the energy issue. Furthermore, we study the hairy black hole solutions 
in which the abelian gauge field breaks the symmetry near the horizon. In the holographic 
picture, it is equivalent to a second order phase transition. Explicitly we show that there exist 
a critical temperature which is a function of the Lifshitz parameter z. The system below the 
critical temperature becomes superconductor, but the critical exponent of the model remains 
the same of the usual holographic superconductors without the higher order gravitational 
corrections, in agreement with Gainsbourg-Landau theories. 
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1 Introduction 



The anti-de Sitter/conformal field theory (AdS/CFT) correspondence [T] has many useful applica- 
tions in condensed matter physics, especially for studying scale-invariant strongly-coupled systems, 
for example, and low temperature systems near quantum criticality [2]. There are many attempts 
to relate the condensed matter problems to their gravitational duals. Since the high-T c super- 
conductors are shown to be in the strong coupling regime, the BCS theory fails and one expects 
that the holographic method could give some insights into the pairing mechanism in the high-T c 
superconductors. From the (d dimensional) field theory point of view, superconductivity is charac- 
terized by condensation of a generally composite charged operator O in low temperatures T < T c . 
In the gravitationally dual (d + 1 dimensional) description of the system, the transition to the 
super conductivity is observed as a classical instability of a black hole in an anti-de Sitter (AdS) 
space against perturbations by a charged scalar field ip. The AdS/CFT correspondence relates the 
quantum dynamics of the boundary operator O to a simple classical dynamics of the bulk scalar 
field ip [3111]. Various holographic superconductors have been studied in Einstein theory [5l|6] or 
extended versions as Gauss-Bonnet (GB) [HIE], Weyl corrected ones [9tiTTj . with magnetic field 
in the bulk action [T2HT8] and even in the non relativistic model of gravity as, for example, in 
Horava-Lifshitz theory [T9||2"0]. In recent years, holographic method have been used to study non- 
relativistic system [21]. In the framework of condensed mater theory, different system show a 
dynamical scaling near fixed points: 

t ->• X z t, Xi^-Xxi, z ^ 1. (1.1) 

As a consequence, instead obeying the conformal scale invariance t — > Xt, Xi — > Xxi , the temporal 
and the spatial coordinates scale anisotropically. The Lifshitz topological black holes and charged 
Lifshitz black holes have previously been discussed in Refs. [22~||2"3] . 

In the present paper we would like to study holographic superconductors in a new background. 
This set up of the s-wave holographic superconductors uses the Lifshitz AdS black hole as the 
gravitational bulk metric in the probe limit. In the Section 2, we introduce a scalar field model 
non minimally coupled with the abelian gauge field in the presence of cosmological constant in 
N dimensions and we derive static, (pseudo-) spherically symmetric (SSS) solutions with various 
topologies. In particular, we will be interested in the AdS-black hole (AdS-BH) Lifshitz solutions. 
In the Section 3, we study the thermodynamical properties of the solutions and obtain the quasi- 
local generalized Misner-Sharp mass as a Killing conserved charge. We also verify the validity of the 
Gibbs equation by using the Kodama-Hayward temperature. In Section 4 we study the hairy black 
hole solutions in which near the horizon the abelian gauge field breaks the symmetry and in Section 
5 we explore the scalar condensation in our Lifshitz black hole solutions by analytical approaches. 
The matching solutions and the critical temperature will be found. Finally, conclusions are given 
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in the last Section. 



2 Bulk asymptotic AdS^ solution 

We will consider the iV-dimensional action of the following model where the scalar field <p is non 
minimally coupled with electromagnetic potential, 

/=/ d N xy/=j[ {R ~ 2 2A) - \d^d^ + V(</>) - ^{F^F, V )\ . (2.1) 

Here, g is the determinant of metric tensor, g^, M. is the space-time manifold, A is a cosmological 
constant, namely A = — (N — 1)(N — 2)/(2L 2 ), L being a length size, and is the electromagnetic 
field strength coupled with scalar field <p as £Exp[A0](i 7Y " / F^ u ) , £ and A being generic constants 
(for example, in the four dimensional Einstein- Maxwell action one has £ = 1/4 and A = 0). Here, 
we use units of fcg = c = h = 1 and denote the gravitational constant k 2 = 87rG^ = 87r(l/M|, z ) 7V 
with the Planck mass of M PL = G^ 1 ^ 2 = 1.2 x 10 19 GeV. The field is also subjected to a potential 

v{4>). 

We look for static, (pseudo-) spherically symmetric (SSS) solutions with various topologies, and 
write the metric element as 

ds 2 = _ e 2 a (r) B{r)dt 2 + __ + ? (2 2) 

where a(r) and B(r) are function of r only and da 2 N _ 2k represents the metric of a topological 
(N — 2)-dimensional surface parametrized by k = 0, ±1, such that the manifold will be either a 
sphere Sat_2 (for k = 1), a torus T^_ 2 (for k = 0) or a compact hyperbolic manifold (for 
k = — 1). In particular, we will be interested in the Lifshitz solutions, where a(r) = logr 2//2 , being 
z the redshift parameter. In this case, for power counting renormalizability in N dimension, if 
we assume that the interaction potential can be expanded as V((p) = £^ =o g m m , where m are 
the polynomial terms of the series and g m suitable coefficients, by the dimensional engineering we 

get m 

[iV+ Z -l-m(iV-l- Z )/2 

[g m \ = [m\ 

Such kind of theory is renormalizable if the couplings have non-negative momentum (here mass) 
dimension, so that we have two possibilities, namely 

2(^-1+^ 

N-l-z 

K = oo, z > N - 1. 

The above constraints are valid for any scalar theory under the Lifshitz scaling of the coordinates. 
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We propose the potential in the following form, 



V(0) = V o e 7 *. 



(2.3) 



where V$ and 7 are generic parameters. Now, by comparing the exponential potential ( 12. 3 h with 
the polynomial form, we observe that in fact K = 00, so that the theory (12.11) is renormalizable 
for z > N - 1. 

With the metric ansatz ( 12.2ft . the scalar curvature reads 



R 



-3B'a' - 2 5a' 2 - B" - 2 £a" — (N — 2) 



l B , + 2 E*_(E^) {k _ B) 

/y ty ry di 



(2.4) 



where the prime index denotes the derivative with respect to r. Where is not necessary, the 
argument of the functions a{r) and B(r) will be dropped. 

Due to the SO(N — 2) symmetry, it is easy to see that the only non vanishing components of 
the electromagnetic field in ^-dimension are 



-a(r) 



Q 



JV-2 



(2.5) 



Q being the electric charge of electromagnetic field multiplied to some suitable dimensional positive 
paramter. By assuming 6 = d(r) and by plugging the above expressions into the action (12. ip . 
making a partial integration, we get the following effective Lagrangian, 



C cS = e a r N ~ 2 



(N - 2)(N - 3)(k - B) (N-2)B' 



2A - ^_ + v** + 2i^$A 



(2.6) 



where, for simplicity, we have putted 6 = 2k, 2 6, Vq = 2k 2 Vq and Q 2 = 2k 2 Q. Making the variation 
with respect to a(r) and B(r) we get the Equations of motion (EOMs), 



B'r - (N - 3)(k - B) 



(N-2) 







(2.7) 



St 

a'-—- - = 0. (2.8) 

This equations correspond to the (t, t) and (r, r) independent components of the usual tensorial 
form of field equations, 



R,u - \g^R + Ag^u = k 2 (T* u + ^(4F va d^ - g^F^F^)) 



(2.9) 



where R^ u is the ^-dimensional Ricci tensor, T^ v is the stress energy tensor of the field 6 and A u 
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represents the electromagnetic potential such that = V M v4„ 



Finally, the equation for <p{r) reads 



a'B + B' 



N -2 



-B 



7W 



2£Ae A *Q 2 

r 2(N-2) 



0. 



As we stated above, we are interested in the following solutions, 



a(r) = log[(r/r ) 2/2 ] 



(2.10) 



(2.11; 



which correspond to the important class of Lifshitz solutions parameterized by red shift z param- 
eter. Here, r$ is a dimensional constant. From Eq. (12. 8ft we get 



(r) = y/z(N - 2) log[r/r ] , 



(2.12) 



so that for renormalizable theory (z > N — 1) we deal with real fields, since in general we take 
N > 3. By using this result, we can solve Eq. (12. 7p as 



B(r) = 2k(N-3) 3 , N+z/2 2r^ 

1 ; z + 2N - 6 (N - 2) (6 - 2A - 2N - z) 



4£Q 2 r 



2^B 



4Ar 2 



(N -2)(Q-2B -2N - z) (N - 2){2N — 2 + z) 
In this equation, C is a free integration constant of the solution and 



(2.13) 



A 
B 



7v / 2 (A-2) + 2, 
Xy/z(N- 2) + 6 - 2JV. 



(2.14) 



-yy/ z(N-2)~ 



^/and r/^^Q 2 ^ 



For simplicity, in the above expression we have also redefined r k —r v &i±u. / 

Q 2 . Note that if we turn off the electromagnetic and field potentials and then take <j) = 0, we recover 
the Reissner Norstrom solution with cosmological constant, 



ds 2 = -B{r)dt 2 -f 
B(r) = k + Cr^ N 



— — + r aa (N _ 2)tk , 



2Ar 2 



+ 



2-3A + A 2 (A-5)(A-2)r 2 



(2.15) 



On the other hand, in the presence of the scalar field, the solution (12 . 1 3[) is acceptable only if the 
Klein Gordon equation (I2.10p for <fi is also satisfied. In this case, the form of the metric is given by 



ds 1 



r J B(r 



2 1 "" +r 2 da 2 

' ' aa N-2,k ■ 



(2.16) 



1 



Let us see for two different cases. 



Absence of cosmological constant. If A = 0, Eq. (12.10j) is satisfied by choosing 

2 



7 : 

A : 



y/z(N-2)' 
(2N - 6) 

y/z(N-2) 

k(N - 3)(N - 2)z + 2Q 2 (2N — 6 + z)£ 

: (2^) 



(2.17) 



The solution becomes 



4[^-3) + 2<^] 
V ; (z-2)(2N -6 + z) 



Cosmological constant A ^ 0. One simple solution is given by 

2N-A 

A 



(2.18) 
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Va 



V<N - 2) ' 
2 

" y/z{N - 2) ' 

2 {n — 2) g 2 e 

_A;(iV-3)(iV-2) 2 g 2 e 
(N-1)Q^-A 



(2.19) 



The solution reads 



B(r) 



2k(N-3) 



2V n 



(2N-6 + Z 

4r 2 (A - Q 2 Q 
~(A^-2)(2A^-2 + z) 



(iV-2) 



+ Cr 3 "^" 2 / 2 



(2.20) 



which asymptotically is a de Sitter/ Anti de Sitter (dS/AdS) solution. Note that if A = 0, 
we get z = —2(N — 2), and, for N > 2, the solution is not accettable, being <fi in Eq. (I2.12p 
immaginary, and the theory becomes non renormalizable. In principle, for any dimension 
N and for any choice of z we can obtain the corresponding Lifshitz solution by setting the 
values of 7, A, Q and Vn appearing in the field lagrangian. 

The solution (12.20D can be asymptotically AdS and will furnish our bachground in studying 
holographic superconductors. 
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3 Black hole solutions and thermodynamics 



The solutions derived in the previous Section may describe charged black holes (BH) in N- 
dimensional manifolds in the presence of scalar field non minimally coupled with electrodynamic 
potential. We recall that event horizon exists as soon as there exists a positive solution r + of 

B(r+) = 0, £'(r+)>0. (3.1) 

We require B(r + ) ^ to avoid the non extremal BHs. 
In the case of solution ( 12. 18ft one has 



Ak(N - 3) + 8Q 2 «e 
C(z-2)(2N -6 + z) 



3-N-z/2 

(3.2) 



and since z > 0, we must require k/C > when N > 3. However, in order to have B'(r + ) > 0, we 
see that only in the topological case k = — 1 we obtain a BH solution. 

Concerning the solutions (I2.15P and ( I2.20p . it is always possible to describe topological BHs by 
making an appropriate choice of the parameters. For example, in the case of solution (I2.20p with 

1 ^_ 3)+ m >0 , ,j r :{j-jy ,<o. 



(2N-6 + z) \ v ' (N-2)J ' (N -2)(2N -2 + z) 

the equation B(r) = has two roots, namely r±, the first one corresponding to the event horizon of 
the black hole (B'(r + ) > 0) and the second one to the (Anti-de Sitter) horizon of the cosmological 
background where the black hole is immersed (B'(rS) < 0). 

In the following, we will assume to deal with solutions whose parameters satisfy the conditions 
(13. ip for some value of r = r + . 

Let us study some physical propriety of these black holes. Since the field equations (12.91) of 
the theory are second order differential equations, we can easily derive a conserved current whose 
charge may be identified with the mass of the BHs. We will follow the approach proposed by 
Wheeler for Lovelock theories [25] in Ref. [26]. For simplicity, we denote with Q^ v the Einstein 
tensor plus the contribute of cosmological constant, namely 

Gp» = R^u - ^g^R + A# Mi , . (3.3) 

By means of time-like Killing vector field K" = (1,0) in iV-dimension, in the case of static metric 
(12. 2 p . one can construct the conserved current 

J^.= G^K U , (3.4) 
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such that 

V v r = 0, (3.5) 
being Q^ v a conserved quantity. A direct evaluation of Jo via ( 13. 3ft leads to 

Jo : = «^ = (^ BW) (^j)^7W] , (3 . 6) 

where 

2A 

W(r) = [t-B(r)]r-'- (Af _ 2)(Af _ ir (3.7) 

The current J M gives rise to a Killing conserved charge. This corresponds to the quasi-local 
generalized Misner- Sharp [27] mass which reads 



B " s(r) s -1? L J dS " = a? y. ""^T - = r w W ' (3 ' 8) 



where £ is a spatial volume at fixed time, dH^ = (dE, 0), and Vjv-2,fe is the N — 2 dimensional 
volume depending on the topology. For example, in the case of the sphere with k — 1, one has 
Vat-2,1 = 2vr( 7V - 1 )/ 2 /r((iV - l)/2), with r(z) the Euler-Gamma function. 

In particular, on shell, that is at the horizon r = r + such that B(r + ) = 0, the quasi local energy 
is identified with the black hole energy E which reads 

(JV - 2)V N ^ / „_ 3 2A 



E := E us{r+ ) = ' ^ - (jv _ 2)(jv _ ir rj • (3-9) 

For example, in the vacuum case of Eq. ( 12.151) with Q = one has 

r N ~ x W(r) = -C y (3.10) 

and the black hole energy reads 



E=- ^ 'Z/^ C, (3.11) 



(N-2)V N ^ ki 
2k 2 

so that the constant of integration is related with the mass of the BH. 

We note that expression ( 13. 9 j) correctly returns the Misner-Sharp mass for asymptotically flat 
solutions (A = 0) in vacuum or in the presence of matter. In particular, for N = 4 and k — 1, 
by explicitly writing the Newton Constant, we get the familiar result E = r+/(2Gn), which 
corresponds, in the vacuum case, to E — —C/{2G^). 

Now, let us show that the Gibbs equation TdS = dE — pdV holds true for the black holes 
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described by the model ( 12. ip . with the Killing energy E obtained below, and the pressure p given 
by electromagnetic and scalar fields. T and S are the temperature and the entropy of the black 
hole, and V is the volume enclosed by the horizon in N — 1 dimensional space. 

For Lovelock gravity the validity of the First Law of black hole thermodynamics has been 
investigated in several places [28H30] . For our static non vacuum case we present a simple derivation 
from the first EOM (12. 7p evaluated on the horizon r = r + , 

(N - 2)V N ^ k B'(r + )r^ (N - 2)V N ^ k ( _ 2A 



2k 2 2k 2 V + (N-2)(N- 1) + 

(3,2) 

2k z r\ 

Here, we have used the fact that B(r + ) = 0. 

All thermodynamical quantities associated with black holes solutions can be computed by 
standard methods. The entropy can be calculated by the Wald method [3TH33] and reads 

2nV N - 2 ,k at_ 2 (q i o\ 

K 

Furthermore, for the static metric (12.21) it is possible to find a characteristic temperature related 
to the event horizon. A natural choice is to take the so called Killing/Hawking temperature [M] 

whose validity may be justified making use of derivations of Hawking radiation [35] or by elim- 
inating the conical singularity in the corresponding Euclidean metric [36] or making use of the 
tunneling method [371[38] . In the above expression, kk denotes the Killing surface gravity, namely 
k k = e Q(r+) J B / (r+)/2, derived from the relation K^V \,K V = k k K u , where K v = (1,0) is the 
time-like Killing vector field. 

However, we would like to remind that in the spherical symmetric, dynamical case, the real 
geometric object which generalizes the Killing vector field is the Kodama field with a related 
conserved current and a related Kodama surface gravity k h [39] . In such natural definition 

of the temperature for dynamical black holes reads as T# = k^/27t, where T# is the Kodama- 
Hayward temperature, in analogy with the static case. 

For spherically symmetric statical metric, the Kodama vector differs from the Killing one as 
K v = e~ a ^K u and as a consequence one finds 

namely Th = e~ a( - r ^Tx and in principle definitions (I3.14p and (I3.15P are different. In vacuum case 
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where a(r) = 0, the two temperatures coincide, but for "dirty" BHs (i.e., in the presence of matter) 
as the ones we are considering, T# ^ Tjj. This is related with the fact that the Killing vector 
cannot be defined unambiguously when the space-time is not asymptotically flat. We stress that 
all derivations of Hawking radiation lead to a semi-classical expression for the black hole radiation 
rate T, 

r = e t k , (3.16) 



in terms of the change AEk of the Killing energy Ex, but if one uses the Kodama energy Eh for 
the emitted particle, one has 

AE H 

(3.17) 



II 



This fact derives by the relationship AEjj = e a ^AEx- From the Eqs. fl3.16p - fl3.17p . one arrives 
at the identity 

^ = ^, (3.18) 

T H T K K ' 

so that the tunneling probability is invariant under different choices of the temperature. 

In Ref. [ID] an attempt to identify the mass of static BHs in modified theories of gravity as 

the integration constant which appears in the vacuum solutions has been done. This result has 

been derived by the EOMs and seems in favor of the Killing temperature with respect to the 

Kodama- Hay ward one, but here, for our non vacuum solutions, Eq. ( 13. 12ft suggests the use of the 

Kodama temperature ( 13.151) . in the attempt to recover the Gibbs relation. In fact, by making use 

of the BH entropy (I3.13p . one can rewrite Eq. (13.121) as 

T H dSw = dE + pdV , (3.19) 

where E is the BH energy (I3.9p . V is the volume enclosed by the horizon in iV — 1-dimensional 
space, V = VN~2,k rN ~ 1 / {N — 1), and p = + Pem is the working term given by scalar (p^) and 
electromagnetic (pem) field pressures, 

p^ = _y(0) = _y o e* 

ISO 2 

Pem = iF^F muv = -^-. (3.20) 

Here, we have reintroduced Vq = Vo/(2k 2 ) and Q 2 = Q 2 /(2k 2 ). In this case, the Gibbs equation 
holds true. 

In conclusion, we prefer to use the Hayward temperature Th (13.151 ) in analogy with therody- 
namic, by starting from the robust defininitions of the neergy as the charge of a conserved current 
and the entropy via Wald method. 
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4 s-wave Holographic superconductors in probe limit 



Our goal in the following Sections is to apply the black hole solution with the Lifshitz scaling which 
it has been obtained before to study the holographic picture of a high temperature superconductor 
via gauge/gravity duality. In brief, by holographic superconductor we mean a condensed matter 
system under second order phase transition whose physical properties can be described by studying 
the dynamics of gauge field on the black hole background in the bulk. It is proved that by direct 
applying the AdS/CFT conjecture one can interpret the asymptotic behavior of the gauge fields 
on the AdS boundary as the expectation values of some physical operators. The expectation 
values of such scalar operators are dual to the super current in the s-wave high temperature type 
II superconductors. Here s-wave refers to the U(l) gauge field. The gauge field can be SU(2) 
and corresponds to the Yang-Mills fields. Such holographic models are called p-wave, because 
super current is a vector and the condensation happens usually for one homogenous component 
of it. To apply the gravitational model to the superconductors, we will modify our gravitational 
model (12. ip by adding a new matter field subjected to some abelian gauge field. By starting from 
the backgound metric of gravitational action previoulsy studied, we will investigate the scalar 
condensation of the new field and we will show that some phase transition occurs. 

We easily see that Eq. (12.201) may describe an (dS/AdS) black hole solution for our non 
minimally coupled gravity model in N dimensional bulk. We are interested in Anti de Sitter 
solutions. If we want to relate our gravitational system with a strongly correlated system in the 
dual quantum theory, we need to describe the dual quantum operators via CFT. The condensed 
matter system dual to our AdS solution can be addressed by holographic superconductors. In 
what follows, we will study the formation of the hairy BHs. The phenomenon is given by a second 
order phase transition and can be described by the holographic methods of the AdS/CFT. 

At first, we note that the solution (12.201) with A ^ may be asymptotically topological AdS 
and it can be considered as the gravitational part of the holographic superconductor in the bulk. 
In fact, we take the gravity bulk as the charged topological black hole with a non zero, negative 
effective cosmological constant 



4(A - Q 2 Q 
(N - 2){2N - 2 + z) 



KfJ = ~ nZn.r o ■ < ■ (4-1) 



It is very interesting that the 17(1) reduced charge Q is combined with cosmological constant in 
the solution, producing the AdS background. It is useful to rewrite the solution (12.201) as 

B " - w^) mN - 3) + + a ^ n + i • (4 - 2) 
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where 



l eff 




2)(2iV-2 +£ ) 



QK - a 



is the effective lenght scale. We note that for N = 4 and z — 0, solution H4.3[) reduces to the usual 
Schwarzschild-AdS form with l e jj = ^/3/A. In fact, when we turn off the electromagnetic field, 
the form of the effective length scale is the same as the one of an AdS uncharged black hole. 

We want to find the second order phase transition in the bulk theory by studying the boundary 
operators. In particular, we want to study the role of the Lifshitz scaling z related to the critical 
temperature T c and the condensation of the dual operators < (D± > (see also Refs. [4TH4"3] ). As 
a starting point, in order to discuss the superconducting phase via holographic picture, we need 
a scalar field ip(r), with mass upper than the Breitenlohner-Freedman (BF) bound [15], and an 
abelian gauge field *4 M , which is non minimally coupled with the scalar field, so that we modified 
the (12.1 p by introducing a new matter Lagrangian in the following form 



C m = —T^T^ - \D^\ 2 - m 2 ?/; 2 , (4.4) 

where is the covariant derivative, = d^ — iqA^, and J 7 ^ is the electromagnetic field strength 
related to the abelian field. As a consequence, the total action results to be 

Itotai = ~\ f d N xV=g£ m + / d N x^9 P~ 2A) - \d»<j>d^ + V{<j>) - & x<t >{F^F, u ) 
Q J Jm i IK I 

(4.5) 

In the above expressions, q plays the role of bulk electric charge. We stress that it is not the 
physical "electric" charge, but it is the dyonic magnetic one. The matter action (I4.4p is different 
from the bulk action. In fact, we add here the non minimally coupled Maxwell field in order to 
break the U(l) symmetry of the abelian field near the BH horizon. Thus, we will work in the so 
called probe limit q — > oo, ignoring the back reaction, in the normal phase ip = 0. In this case, 
the gravity sector decouples from the abelian one and the background metric can be derived from 
Eq. fl2~T6l) and Eq (12^201) . 

In the probe limit the EOMs for ip(r) and = ^(r)^, 5^ being the Kroenecker delta function 
and 0(r) a general function of r, read in the following forms 

D^D^ip - m 2 ij = 0, (4.6) 
V\7> = iqWD v ^ - $DW\ . (4.7) 

Here, = d^A u — d u A^ and we take ip* = i/j, motivated by the fact that we are free to choose 
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our gauge. The above equations read in terms of the background metric as 



N-2- 



2 I JJ 



2„/,2 



2q z tp 
~B 







i= + ( + — W 



m 2 q 2 r z d 2 \ 



(4.8) 
(4.9) 



For N = A and z = the model is perfectly described by numerical methods [46] . We are interesting 
in the cases of z ^ and iV > 3. 

We need to introduce in Eqs. (14 .81) - (14 .9p the horizon radius r + . We rewrite the solution in the 
following form, 



B(r) = k 



. 3-iV-f 



r+ y 



, 3-JV-f 



(4.10) 



where 



2N 



( 2k(N - 3) + -J^ 



(4.11) 



It is more convenient to work in terms of the dimensionless parameter y{r) = r + /r, such that 
y(r+) = 1 and at the infinity y(r — > +oo) — > + @ In this case, the equations of motion (14.81) and 
(14. 9 p can be expressed as: 



+ N-4\ , 2r^ 2 



y \ B 



B 2 



if) = 



(4.12) 
(4.13) 



and the metric function reads 

B( 3 j) = k(l-y N -^ 



5 + 



eff' 



) \y -y 



(4.14) 



Now, the prime denotes the derivative with respect to y — y(r). We have obtained the basic set 
up for the holographic superconductors. 



5 In the literature usually the symbol is z, but here we kept z for Lifshitz scaling and we introduced y as the 
dimensionless radial coordinate. 



12 



5 Critical temperature and condensation values by match- 
ing method 

According to the AdS/CFT correspondence, we can interpret (O), where O is the operator dual 
to the scalar field, as the expectation value of the dual current operator on the AdS boundary 
r — > +00. Thus, we are going to calculate the condensate (O) for fixed charge density. 
Regularity at the horizon, namely y = 1, requires 

0(1) = 0, f(l)B'(l)=rX^(l). (5.1) 

We want to find approximate solutions around the horizon and asymptotically AdS limit, 
y = 1 and y = 0, using Taylor's expansion, then we want to connect these solutions in an arbitrary 
matching point y m between y = 1 and y = 0. In principle, we may do the computations by 
numerical methods using shooting algorithm, but in this paper we would like to keep the level of 
the analytical approach. At first, we calculate B'{1) (= (dB(r + ) / dy) directly from Eq. (I4.14p . 



B>(1) = -A (J V - 3 + |) + + 1 - |) - (5.2) 



It is easy to see that Eq. (15 .2p is related to the Kodama temperature (I3.15P by B'{1) = —4irr + TH- 
We will use (15. 2 p to construct the series solutions for our topological holographic superconductor 
in the probe limit. 

5.1 Solution near the BH horizon 

We expand and ip near the black hole horizon at y = 1 as 

1 

2 
1 



<f>{y) = 0(l)-0'(l)(l-y) + -0"(l)(l-y) 2 + ---, (5.3) 



tf>(y) = ^(l)-^'(l)(l-y) + ^"(l)(l-y) 2 + ---. (5.4) 

From the boundary condition, we know that (f)(1) = and for simplicity we put a := — 0'(1) < 
and b := ^(1) > for the positivity of (j)(y) and ip(y). 

In order to discuss phase transition near the critical points, we need just to keep the second 
order terms in those series. 

First, we compute the 2 nd order coefficient of <fi by using Eq. (14.121) , 



mi = -<b'(i) 



,z \ 2r 2 ,b 2 i , 
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In this case Eq. ( 15.31) reads 



1 



M) = a(l-y) + ^ -(| + JV-4) + ^](l-y) s 



We can calculate the 2nd derivative of ijj from ( 14.131) in the same way, 



ba 2 (r + ) 2 " z 



and write the following series solution for scalar field ip ( 15. 4 p 



i i 
m r± 



a 2 r 2 r z 



1 B'(l)^ ^ 2(5'(1)) 2 
In the above equations, B'{1) is given by Eq. ( 15 .2[) . 



(5.6) 



(5.7) 



(5i 



5.2 Solution near the asymptotic AdS region 

The asymptotic regime of the metric function B(y) in the AdS boundary is completely independent 
on the dimension of the spacetime TV. It is trivial to recover the following asymptotic behavior of 
the metric 



B(y) ~ y 



(5.9) 



Here we take z > and N > 3 to avoid the problems of convergence in AdS /C FT. Consider now 
the weak field behavior (r — > +oo) of which depends on the value of dimension N. It is easy to 
show that this behaviour completely changes at the critical dimension iV c = 3 — z/2, namely 



<i>{y) 



Cl 



N -3 



y 



;+N-3 



+ c 2 , N^N C 



Cl log(y) + c 2 , A^ = A^ C . 



(5.10) 
(5.11) 



We are interested in the fields with fall off behaviors near y = 0, so that we take N ^ N c case as 
the physically acceptable solution. It is useful to write (I5.10p in terms of the radial coordinate r, 



ciri 



|+iV-3 



+ N -3)r 



\+N-3 



+ C 2 . 



(5.12) 



By writing this solution in terms of the dual physical quantities chemical potential \x and charge 
density p, we obtain 



(5.13) 
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where 



P 



i+N-3 



+ N-3) 



P 

|+iV-3 



(5.14) 



The second condition derived from the fact that 4>{r + ) = 0. Also for the scalar field, near the AdS 
boundary we can write 



<0,> < CL > 



(5.15) 



where < + > and < CL > are the operators on the boundary. Furthermore, to have a normal- 
izable solution, we will set the source < CL >= 0. By rewriting it in terms of y with the fall off 
behavior just for operator + we get 



— A—y ■ 



(5.16) 



Here 



■\{N - 1) ± yJ(N - l) 2 + Am% f ] 



(5.17) 



5.3 Matching and phase transition 



In this Section, we will connect the solutions ( 15 .6p and (15. 8p with (I5.13P and ( I5.16P at some 
matching point y = y m . In order to connect those solutions smoothly, we require the following 
four conditions: 



/' 



P z/2+N-3 

■ y^ u 



r 



,z s 2r\b 2 

^2 ' B' 1 



p(z/2 + N-3) ^ z/2+N 



r 

<o + > 



z/2+N-3 ^ m 



- 1 



,z . 2rlb 2 i / 



-v, 



r , 



1 - 



2 2 



Ml] 



[i -yn 



A+<0+> 



yt~ l = b 



i i 

m rl 



+ 



a 2 r 2 ~ z 



B'(l) (S'(l)) 



a 2 r 2 + ~ z 
2(B'(1)) 2 

i - y m 



(5.18) 
(5.19) 
(5.20) 
(5.21) 



By combining Eqs. ( I5.18p -( l5.19p we can eliminate ab 2 and one has 



2 p ( (f + AT - 3) 2 /_f+ JV - 3 - (-5 + f + iV)y m f +7V " 2 )r + -f- w + 3 + 2(1- 2/m ) ?/m a 



4y„ 



V2 
2r + ^ 



PZ/ m i +7V - 3 (| + N-3)-((l+N- A)y 



N + 5)y m ar 



l+N-3 



B'if 



\+N-3 



?/m(l - yn 



(5.22) 



(5.23) 
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The above relations allude to the phase transition, namely, given p, fi has a maximum value when 
we assume the non-trivial solution 6^0. Now we can reveal the phase transition in our simple 
system. In order to evaluate the expectation value of the operator < + >, we eliminate the a 2 b 
term from (15.201) and f 1 5 . 2 1 1) and obtain 



r + A b(m 2 r + 2 (y m - 1) + 2)y, 



l-A 



((A - 2)y m - A) 

For non- vanishing b, we can compute < + > from Eqs. (I5.20I) - (I5.21I) and one gets 



(5.24) 



a = \B'(V 



2A 2 -m 2 r 2 + (y m + 2A(l-y n 



^r 2 + z [A(l-y m ) 2 + 2y m (l- ym ) 
By plugging this result in Eq. (I5.23P we derive 



(5.25) 



2r A 



V2B'(1) 



\ (i -ym)Sr. 



'i+N-3 



+ N- 3)pyl +N 4 + y/2: 



S = \B'(V 



(m 2 r + 2 (A - l)y m - A (-1 + m 2 r + 2 )) 
r + 2 -(y m -l)((A-2)y m -A) 



(5.26) 
(5.27) 



By using the density (I5.14p . one has that (0+) can be expressed as 

(0 + ) = 

where the new function T is defined as 



y m l -^{r + 2 {-l + y m )m 2 + 2 )r + A r= 
((A - 2)y m - A)r + 



(5.28) 



r 



( - v^((f + N- A)y m - f - N + 5)S r + § +Af - 3 + (§ + #- % m § +7V - 4 p) ■ 



Er + * + "- s (l-&») 
Now we can write T in the following equivalent form 

Tr-T H 



(5.29) 



(5.30) 



H 
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where T# is the Kodama temperature (13. 15ft and A is a function of {r + ,A,y m }. The critical 
temperature T c is defined as 



f+TV-3 



(f+jV-3)jfa 

VP 4V2y m vrl +N - 2 (-(l + N- A)y m + § + TV - 5) 



r + 2 -(l-?/ m )((A-2)?/ m -A 
\ — m 2 r + 2 (A — l)y m — A (1 — m 2 r 



(5.31) 



Note that if we choose y m = 1/2, z = and Af = 4, the critical temperature behaves like T c ~ ^/p, 
according to the literature about s-wave holographic superconductors. This behavior will change 
only in the presence of the higher order corrected backgrounds like Weyl's models for s-wave. For 
example for Gauss-Bonnet and Weyl corrections it reads as T c = ^fp [H]. 
Finally, using Eq. (15.311) we have 



(0 + ) - ^:J H • (5-32) 

V J-H 

The last expression has been written near the critical point T c . 

We see that (0 + ) is zero at T# = T c , which is the critical point, and condensation occurs at 
T H < T c . The behaviour of (0 + ) oc (1— Th/Tc) 1 ' 2 is also recovered, in agreement with the literature 
about the argument. The system below this critical temperature becomes superconductor, but the 
critical exponent of the model remains the same of the usual holographic superconductors without 
the higher order gravitational corrections. 



6 Discussions 



In this paper we have considered a holographic model for a non-relativistic system showing su- 
perconductivity. We have used a black hole background which comes from a scalar field model 
non minimally coupled with the Abelian gauge field in the presence of cosmological constant A 
in N dimensions, and we have studied analytically holographic superconductors in this new kind 
of asymptotic AdS solutions. We have considered static, (pseudo-) spherically symmetric (SSS) 
solutions with various topologies in two different cases, A = and A ^ 0. We have obtained the 
quasi-local generalized Misner-Sharp mass as a Killing conserved charge. This quasi-local energy 
at the horizon r = r + is identified with black hole energy. Then we have derived the Wlad entropy, 
Killing-Hawking temperature and Kodama-Hayward temperature of black hole solutions. These 
temperatures are in principle different. We have shown that for our non- vacuum solution, the first 
law of black hole thermodynamics holds true by making use of Kodama temperature. After that 
we have studied the hairy black hole solutions in which near the horizon the Abelian gauge field 
breaks the symmetry. In the holographic picture, this symmetry breaking is equivalent to a second 
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order phase transition near the horizon. We also have analytically solved the system in the probe 
limits, near horizon and asymptotic region. We have found that there is also a critical temperature 
which is a function of the Lifshitz parameter z and under such temperature a condensation field 
appears. 

Acknowledgments 

We would like to thank Dr. Yanyan Bu and Prof. R. B. Mann for advices. L.S. thanks also Prof. 
S. Zerbini for valuable suggestions. 

References 

[I] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998). 

[2] C. P. Herzog, J. Phys. A 42 (2009) 343001 ; S. A. HartnoU, Class. Quant. Grav. 26 (2009) 
224002. 

[3] S.S. Gubser, I. R. Klebanov and A. M. Polyakov, Phys. Lett. B 428 (1998) 105. 

[4] E. Witten, Adv. Theor. Math. Phys. 2 (1998) 253. 

[5] S. A. HartnoU, C. P. Herzog and G. T. Horowitz, JHEP 12, 015 (2008). 

[6] G. T. Horowitz and M. M. Roberts, Phys. Rev. D 78, 126008 (2008). 

[7] R. Gregory, S. Kanno and J. Soda, JHEP 10, 010 (2009). 

[8] Y. Brihaye and B. Hartmann, Phys. Rev. D 81, 126008 (2010). 

[9] D. Momeni, M.R. Setare JarXiv:1106.0431l 

[10] D. Momeni, N. Majd, R. Myrzakulov JarXiv: 1204. 12461 . 

[II] D. Momeni, M. R. Setare, R. Mvrzakulovl arXiv: 1209.3~l0l 
[12] M. R. Setare, D. Momeni JarXiv:1106.1025l 

[13] D. Momeni, E. Nakano, M. R. Setare, W.-Y. Wen, larXiv: 1108.43401 
[14] M. R. Setare, D. Momeni, R. Myrzakulov, M. Raza JarXiv:1210.1062l 

[15] O. Domenech, M. Montull, A. Pomarol, A. Salvio, P. J. Silva,JHEP 
1008:033.f2010) JarXiv:1005.1776[ 



18 



[16] M. Montull, O. Pujolas, A. Salvio, P. J. Silva, larXiv: 1105.53921 
[17] M. Montull, O. Pujolas, A. Salvio, P. J. Silva, larXiv: 1202.00061 
[18] A. Salvio, larXiv:1207.3800l 

[19] R.-G. Cai and H.-Q. Zhang, Phys. Rev. D 81, 066003 (2010). 

[20] D. Momeni, M. R. Setare and N. Majd, JHEP 1105, 118 (2011). 

[21] S. Kachru, X. Liu and M. Mulligan, Phys. Rev. D 78, 106005 (2008). 

[22] R.B. Mann, JHEP 0906 (2009) 075 , larXiv:0905.lT36l 

[23] M. H. Dehghani, R. B. Mann and R. Pourhasan , larXiv:1102.0578l 

[24] M. Visser, Phys.Rev.D80:025011(2009)larXiv^09m0590J 

[25] D. Lovelock, J. Math. Phys. 12 , 498-501 (1971). 

[26] J. T. Wheeler, Nucl. Phys. B 268, 737 (1986). 

[27] C. W. Misner and D. H. Sharp Phys. Rev. 136 B571 (1964). 

[28] R. C. Myers, J. Z. Simon, Phys. Rev. D 38, 2434 (1988). 

[29] A. Paranjape, S. Sarkar and T. Padmanabhan, Phys. Rev. D 74, 104015 (2006). 

[30] H. Maeda, M. Nozawa, Phys. Rev. D 77, 064031 (2008) |arXiv:0709. 11991 [hep-th]]; 

H. Maeda, S. Willison and S. Ray, Class. Quant. Grav. 28, 165005 (2011) (arXiv: 1103.41841 
[gr-qc]]. 

[31] R. M. Wald, Phys. Rev. D 48, 3427 (1993). 
[32] M. Visser, Phys. Rev. D 48, 5697 (1993). 

V. Faraoni, Entropy 12, 1246 (2010) jarXiv: 1005.23271 [gr-qc]]. 



[34] S. W. Hawking, Nature 248 30 (1974); 
Commun. Math. Phys. 43 199-220 (1975). 

[35] M. Visser, Int. J. Mod. Phys. D 12, 649 (2003) |hep-th/0106TIT] . 

[36] R. Di Criscienzo, M. Nadalini, L. Vanzo, S. Zerbini and G. Zoccatelli, Phys. Lett. B 657 
(2007) 107 [arXiv:0707.4425l [hep-th]]; 

R. Di Criscienzo and L. Vanzo, Europhys. Lett. 82, 60001 (2008) |arXiv:0803.0"435l [hep-th]]; 
S. A. Hayward, R. Di Criscienzo, L. Vanzo, M. Nadalini and S. Zerbini, Class. Quant. Grav. 26, 
062001 (2009) [arXiv:0806.00T4l [gr-qc]]; Class. Quant. Grav. 27 015006 (2010). 

19 



[37] M. K. Parikh and F. Wilczek, Phys. Rev. Lett. 85, 5042 (2000). 

[38] M. Angheben, M. Nadalini, L. Vanzo and S. Zerbini, JHEP 0505, 014 (2005) 
|arXiv:hep-th/0503081] ; 

M. Nadalini, L. Vanzo and S. Zerbini, J. Phys. A 39, 6601 (2006) |arXiv:hep-th/0511250] . 
[39] S. A. Hayward, Class. Quant. Grav. 15, 3147 (1998). 

[40] G. Cognola, O. Gorbunova, L. Sebastiani and S. Zerbini, Phys. Rev. D 84 (2011) 023515 
jarXiv: 1104.28141 [gr-qc]]. 

[41] E. J. Brynjolfsson,U. H. Danielsson, L. Thorlacius,T. Zingg. arXiv:0908.26111 

[42] E. J. Brynjolfsson, U. H. Danielsson, L. Thorlacius , T. Zingg . larXiv: 1004. 5566, 

[43] Y. Bu, Phys. Rev. D 86, 046007 (2012). 

[44] R. Gregory, S. Kanno, J. Soda JarXiv:0907.3203l . 

[45] P. Breitenlohner , D. Z. Freedman, Annals Physics 144, 249 (1982). 

[46] S. A. HartnoU, C. P. Herzog, G. T. Horowitz, Phys. Rev. Lett. 101, 031601, (2008). 

[47] X.-M. Kuang, W.-J. Li, Y. Ling , larXiv: 1008.40661 . 



20 



